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ABSTRACT 


The  prediction  of  wave  characteristics  generated  by  variable 
winds  under  fetch-and  duration-limited  conditions  is  considered 
by  using  a  parametric  approach  and  a  simplified  form  of  the 
radiative  transport  equation.  The  wind-generated  sea  surface  is 
characterized  in  terms  of  an  equilibrium  range  spectrum  with  a 
well-defined  low  wave  number  cut-off.  By  relating  the  local 
wind  velocity  to  the  rate  of  increase  in  wave  energy  and  incor¬ 
porating  the  refractive  and  frictional  effects  of  transitional 
water  depths  into  the  formulation,  governing  equations  are  derived 
to  predict  the  development  of  a  wind-wave  field  in  terms  of 
integral  properties  of  the  spectrum  such  as  mean  period,  direction, 
and  significant  wave  height.  Under  deep  water  conditions, 
governing  equations  admit  simple  closed  form  solutions  which  are 
shown  to  compare  very  favorably  with  observational  data  as  well 
as  various  wind-wave  prediction  formulae  currently  available. 
Consequently,  they  are  numerically  solved  first  for  the  case  of 
wind-waves  in  uniform  water  depth  and  then  for  a  circular  basin 
with  variable  bathymetry  as  an  illustrative  general  case. 
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INTRODUCTION 


The  field  of  wind-generated  waves  has  been  the  subject  of 
intensive  research  effort  over  the  last  three  decades  and,  as 
a  natural  consequence,  a  great  deal  of  progress  has  been  made 
on  various  aspects  relevant  to  understanding  the  physics  of 
wind  waves,  their  generation,  propagation  and  dissipation.  An 
excellent  review  of  recent  advances  in  the  theory  of  wind  waves, 
their  observation,  measurement  and  prediction  is  given  by 
Barnett  and  Kenyon  (2).  It  appears  that,  in  spite  of  these 
advances  achieved,  the  present  state  of  knowledge  is  still  in 
need  of  considerable  improvement,  and  various  gaps  in  theory, 
observation  and  instrumentation  remain  to  be  filled  through 
future  research. 

One  ot  the  primary  applications  of  wind  wave  research  is 
wave  forecasting.  The  complete  problem  is  a  field  radiation 
problem  in  which  the  directional  energy  spectrum  of  a  wave  field 
is  governed  spatially  and  temporally  by  an  energy  balance  or, 
radiative  transfer  equation  (1,2,3,10,11,21,22,23).  The  solution 
of  the  radiative  transfer  equation  requires  integration  ot  various 
non-linear  source  terms  representing  the  action  of  local  winds, 
wave-wave  interactions,  dissipative  forces,  etc.  along  wave 
rays  which  are  the  propagation  paths  given  by  the  Hamilton- 
Jacobian  equations.  However,  the  integration  is  not  straight¬ 
forward,  since  the  source  terms  are  non-linear.  This  means  that 
different  rays  are  coupled,  making  the  solution  of  the  problem 
quite  difficult  and  impractical  under  general  conditions.  As  a 
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result,  most  of  the  explicit  applications  such  as  those  demon¬ 
strated  by  Hasselmann  and  Collins  (11),  Resio  and  Vincent  (22, 

23),  and  others  (1,2,3,7,14,19)  are  either  based  on  simplified 
approximations  of  various  source  terms  or  consider  deep  water 
conditions.  At  present,  there  is  a  large  number  of  such 
computer-based  discrete  spectral  models  for  forecasting  waves 
(see,  e.g.,  Refs,  7,14,19). 

A  second  approach  which  avoids  the  intricate  complexity 
of  discrete  spectral  models  requiring  large  computers  and  com¬ 
plex  programming  is  the  parametric  specification  and  fore¬ 
casting  of  wind  waves.  This  approach  is  based  on  one¬ 
dimensional  generation  models  which  parameterize  some  integral 
property  of  the  wave  spectrum  such  as  wave  height,  or  wave 
period  in  terms  of  wind  speed,  fetch  and  wind  duration.  In¬ 
cluded  in  this  group  are  the  original  Sverdrup  and  Munk  model 
(28),  and  its  subsequent  versions  improved  by  Bretschneider  (4, 

25),  and  Ijima  and  Tang  (12,25)  based  on  quasi-empirical ,  quasi- 
theoretical  arguments,  and  the  methods  developed  by  Pierson, 

Neuman  and  James  (20)  and  Wilson  (30).  More  recently,  Hasselmann 
et  al.  (10)  have  proposed  a  parametric  model  based  on  the 
assumption  of  the  invariant  spectral  shape  of  the  JONSWAP  spectrum. 
The  validity  of  this  model  depends  on  the  dominance  of  the 
self-stabilizing  effects  of  wave-wave  interactions  as  compared 
to  external  sources  of  energy.  A  comparison  of  various  para¬ 
metric  models  is  given  by  Resio  and  Vincent  (22). 

Parametric  and  discrete  spectral  models  have  been  experi¬ 
mented  with  mostly  under  open  ocean  conditions.  Except  for  the 
recent  applications  of  discrete  spectral  models  to  the  Great  Lakes 
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by  Kesio  and  Vincent  (22,23) ,  and  the  parametric  procedures 
developed  by  Bretshneider  (4,25),  and  Ijima  and  Tang  (12,25) 
tor  conditions  of  unlimited  wind  duration  in  uniform  water 
depth,  there  has  been  a  lack  of  effort  to  extend  either  ot 
these  approaches  to  wave  field  predictions  involving  variable 
winds  in  areas  bounded  by  irregular  shorelines.  Various  pra¬ 
gmatic  preocedures  and  definitions  of  effective  fetch  and 
effective  duration  have  been  used  to  incorporate  such  effects 
into  the  nearshore  applications  (see  e.g.,  Ref.  24). 

The  principal  motivation  here  is  to  develop  a  method  ot 
forecasting  nearshore  wave  characteristics  generated  by  variable 
winds  under  tetch-and  duration-limited  conditions  by  using  a 
parametric  approach  and  a  simplified  form  of  the  radiative 
transfer  eqaation.  It  is  assumed  that  the  wind  sea  is  characte¬ 
rized  in  terms  of  an  equilibrium  range  spectrum  with  a  well- 
defined  low  wave  number  cut-off.  By  relating  the  local  wind 
velocity  to  the  rate  of  increase  in  wave  energy  and  incorporating 
the  refractive  and  frictional  effects  of  transitional  water 
depths  into  the  formulation,  governing  equations  are  derived  to 
predict  the  development  of  a  wind-wave  field  in  terms  of  the 
integral  properties  of  the  spectrum  such  as  mean  period  and 
significant  wave  height.  Various  incompletely  understood  or 
complex  mechanisms  such  as  wave-wave  interactions,  turbulent 
dissipation  due  to  wave  breaking,  wave-current  interactions, 
percolation,  etc.  are  neglected.  Under  deep  water  conditions, 
governing  equations  admit  simple  closed  form  solutions  which 


are  shown  to  compare  very  tavorably  with  observational  data  as 
well  as  various  wind-wave  prediction  formulae  currently  avail¬ 
able.  Consequently,  they  are  numerically  solved  first  for  the 
case  of  wind-waves  in  uniform  water  depth  and  then  for  a  circular 
basin  with  variable  bathymetry  as  an  illustrative  general  case. 

ASSUMPTIONS  AND  FORMULATION 

With  the  neglect  of  current-wave  interactions,  the  energy 
balance  of  wind  generated  waves  can  be  expressed  as  (see, 

Kef.  21,  page  147) 

d¥/dt  =  (3/3 1  +$.$)¥=  T  (1) 

where  ¥(k)  ■  directional  spectral  density  over  a  two-dimensional 
wave  number  space  &  =  (kj,k2);  $  =  group  velocity  associated 
with  the  spectral  component  with  wave  number  vector  ic;  $  = 
horizontal  gradient  operator;  and,  r  represents  various  source 
terms  from  any  process  that  transfers  energy  into  and  out  of 
the  spectral  component  ¥(k).  For  a  detailed  discussion  of  the 
general  form  of  r,  the  reader  is  referred  to  Refs  10  and  23. 

Eq.  (1)  implies  that  the  rate  of  change  of  spectral  density 
along  the  propagation  path  of  a  wave  group  is  T,  which  is  in 
general  a  complex  non-linear  function  of  the  local  wind  and  the 
spectral  density  ¥  itself.  The  propagation  path  is  given  by  the 
Hamilton-Jacobian  equations 

dx£/dt  =  3o/3ki  =  Vi  (2) 
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(3) 


dk^/dt  =  -3o/3x£ 

where  (i=l,2)  represents  a  horizontal  coordinate  system  in 
llu*  still  water  plane; 

o  =  (gk  tanh  kd)^  (A) 

is  the  radian  frequency  with  g  =  gravitational  constant;  and, 
d  =  local  still  water  depth.  Eqs.  2  and  3  describe,  respectively, 
the  wave  ray  associated  with  the  spectral  component  ¥(&)  and  the 
variation  of  &  along  that  ray.  Therefore,  assuming  that  the 
functional  form  of  r  is  known,  the  solution  of  the  wind-wave 
problem  involves  solving  first  the  refraction  problem,  i.e., 

Eqs.  2  and  3.  Once  the  rays  are  determined,  Eq.  1  is  integrated 
along  each  ray.  As  previously  mentioned  in  the  introduction, 
the  integration  is  not  straightforward,  since  I-  is  non-linear 
and  wave  rays  arc  coupled.  Consequently,  in  order  to  overcome 
some  of  these  difficulties,  we  will  appeal  to  various  simplifying 
approximations  as  follows. 

Functional  Form  of  ¥(it).  -  For  wind  generated  waves,  the 
spectral  density  function  V  is  observed  to  have  an  equilibrium 
range  with  a  well-defined  low-wave-number  cut-off  (13,21,29). 
Combining  the  dimensional  reasonings  of  Phillips  (21), 
Kitaigorodskii  et  al.  (13),  and  Thornton  (29),  it  can  be  shown  that 

Y(ic)  =  (3/2ir)  n  k”^  ;  k  £  km,  |a|  <  ir 
—2 

where  S  =  1.31  x  10  ;  o  =  direction  of  the  saturated  components; 

k  =  low-wave-number  cut-off;  and, 
m 
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n  =  %  {1  +  (2kd/sinh  2kd)} 


(5) 


The  spectral  density  as  a  Junction  ot  wave  number  moduli  or, 
the  distribution  of  the  mean  square  surface  displacement  n2 , 
regardless  of  direction,  is  given  by 

F(k)  = /  V(S)  k  da  -  Snk”3  ;  k  >  km  (6) 

-IT 

If  it  is  assumed  that  the  spectral  density  ¥  is  mostly  con¬ 
centrated  in  the  direction  of  the  spectral  peak  wave  number 
l<m,  F(k)  is  approximately  equal  to  the  integral  of  4* (k)  along 
lines  perpendicular  to  the  direction  of  lcm  (see,  e.g..  Ref.  21, 
page  130) . 

One-Dimensional  Form  of  Energy  Balance.  -  Considering  the 
direction  ot  the  spectral  peak  wave  number  as  a  reference, 
a  contracted  form  of  the  energy  balance  equation  can  be  obtained 
by  integrating  Eq.  1  along  lines  perpendicular  to  icm.  Provided 
that  ¥(&)  is  mostly  concentrated  around  icm,  the  contracted  form 
can  be  approximated  as 

dF/dt  =  (3/3 t  +  $  .  $)  F(k)  =  S  (7) 

Functional  Form  of  Source  Term  S.  -  Among  various  possible 
processes  that  contribute  to  the  source  term,  most  of  the  dis¬ 
cussion  will  be  centered  on  the  generative  action  of  overwater 
winds  and  the  dissipative  effects  of  turbulent  bottom  friction. 
Guided  by  previous  studies  by  Miles  (17),  Snyder  and  Cox  (26), 


Man ton  (15,16),  Barnett  (1),  Barnett  and  Wilkerson  (3),  and 
Hasselmann  and  Collins  (11),  the  functional  form  of  S  is  taken 


S  =  F(km)  (  ft  -  *£) 


where 


ft  -  Y  (Pa/Pw)  <km  •  W  -  om) 


with  (Pa/Pw)  =  ratio  of  air  to  water  density;  i5  -  wind  velocity 
at  10  m  above  still  water;  and  y  -  a  constant  (<1)  representing 
the  fraction  of  wind  input  energy  that  is  transferred  down 
the  spectrum  to  build  up  the  lower  wave  number  end  of  the 
spectrum  F(k).  Hence,  it  is  assumed  that  an  equilibrium  exists 
such  that  wind  input  balances  the  energy  losses  due  to  breaking 
away  from  the  spectral  peak  towards  higher  wave  numbers.  The 
function  <|>f  represents  the  dissipative  effect  of  turbulent 
bottom  friction.  For  a  Gaussian  sea  state  in  which  the  spectrum 
is  narrow-beam  centered  around  i^,  the  general  form  of  this  term 
(see,  e.g.,  Refs.  11  and  23)  can  be  replaced  with  the  unidirec¬ 
tional  approximation 


♦f  “  2  Cf  (g  k2/o2  cosh2  kmd)  <|u|> 


where  Cf  =  an  empirical  drag  or  friction  coefficient;  and, 
“  * 


< | u | >  »  { (2/ir)  /  (gk/o  cosh  kd)2  F(k)  dk}* 


7 


GOVERNING  EQUATIONS 

Substituting  Eq.  6  into  the  left  hand  side  of  Eq.  7, 
we  obtain 

dF/dt  =  6nk^3  {-3k”1  dkm/dt  +  n_1(3n/3km) 

(dkm/dt)  +  n-1 (3n/3d)  (dd/dt)}  (12) 

where 

3n/3k  =  d(l-2k  d  ctnh  2k  d)  /  sinh  2k  d  (13) 

m  m  m  m 

3n/3d  =  k  (l-2k  d  ctnh  2k  d)  /  sinh  2k  d  (14) 

m'  m  m  m 

Similarly,  substitutions  of  Eqs.  4  and  6  into  Eqs.  6,10,  and  11 
yield 


a  =  Y(pa/pw>  kro  w0  {w  cos<6-8u)  -  (on/knl  w0))  (15) 

£  =  (8Cj  k^/sinh  2kmd)  {^gdir”1  /  (sinh  v  +  v) 

2k  d 
m 

(v  sinh  v)“2  d v}'“  (16) 

where  WQ  -  a  reference  10  m  -  wind  speed  at  a  fixed  x  and  t; 
w  =  W/WQ,  a  dimensionless  function  representing  spatial  and 
temporal  variation  of  wind  field;  and,  0»6W  -  directions  of 

and  ft,  respectively,  relative  to  a  fixed  reference  axis. 


We  can  now  introduce  the  dimensionless  variables 


K  *  “mWo/8 

1  *  «mv*  <l7> 

t  =  gt/w0 

x  =  g  x/W^ 

d  =  gd/W l 
so  that 

C  =  Cro/W0  =  ojn/kinW0  =  (tanh  k  d/k)*  (18) 

and 

T  =  g  Tm/2sW0  =  (k  tanh  kd)'^  (19) 

represent,  respectively,  the  dimensionless  wave  speed  and 
period  associated  with  the  spectral  peak.  Likewise,  the  dimen¬ 
sionless  significant  wave  height  can  be  shown  to  be 

y*  00 

H  =  «Hs/Wo  =  4  <  /  F(v)  dv  )%  (20) 

J  k 

Finally,  by  substituting  Eqs.  12  through  17  into  Eq.  7  and 
solving  for  dk/dt,  we  obtain 
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(21) 


dk/dt  =  (a/at  +  £  •  i)  k  *  (sx  -  s2  -  s3)  s"1 

where  $  =  dimensionless  horizontal  gradient  operator; 

5  =  n  C  jc  (22) 

with  ic  =  jc/k  *  (cos  0,  sin  0)  ;  and, 

S1  =  Y^pa^pw^  n  -2  *w  cos(0”ew)  “  (tanh  £  d/k)^} 


sinh  2  k  d  (23) 

00 

52  =  8  Cf  n  k^  (3d  j  (sinh  v  +  v) 

J2k  d 

(v  sinh  v)"2  dv}"5  (24) 

53  =  k2  (1  -  2  k  d  ctnh  2  k  d)  d(d)/dt  (25) 

=  kd  (1-2  k  d  ctnh  2  k  d)  -  3n  sinh  2  k  d  (26) 

Hence,  the  variation  of  the  modulus  k  is  described  by  Eq.  21 
along  the  wave  ray  given  by 

dx/dt  *  f  --  n  C  k  (27) 

What  remains,  therefore,  is  the  derivation  of  an  expression 


A 

which  will  describe  the  variation  of  k  or  6  along  the  same  ray 


due  to  refraction.  This  easily  follows  if  we  note  from  Eqs.  3,4 
and  17  that 

dk/dt  =  d (k  k)/dt  =  -  (o  k/sinh  2kd)  ?d  (28) 


Solving  this  for  dk/dt  and  taking  the  dot  product  of  the  re- 

A 

suiting  expression  with  m  =  (-sin  9,  cos  0)  immediately  gives 

de/dt  =  -(o/sinh  2  k  d)  m  ,  V  d  (29) 

In  summary,  the  complete  problem  requires  the  solution  of 
the  system  of  Eqs.  21,27,  and  29  for  a  specified  wind  field, 
i.e.,  WQ,  w,  and  6W,  with  appropriate  initial  and  boundary  con¬ 
ditions  on  £.  Then,  the  estimation  of  wind-wave  parameters 
C,  T  and  H  follows  immediately  from  Eqs.  18,19,  and  20,  respectively. 

COMPARISON  WITH  EXISTING  THEORIES  IN  DEEP  WATER 

Under  deep  water  conditions  with  constant  wind  direction 
@w,  6  =  ©w,  and  no  refraction,  it  can  be  verified  that  Eqs.  21 
and  27  become,  respectively 


dC/dt  =  yQ  C"1  [  w(x ,  t)  -  c] 

(30) 

dx/dt  *  h  C 

(31) 

where  y  *  (1/6)  y  (p  / p  )  ;  and  use  has  been  made  of  k  * 

U  8  W  ™ 

—2 

C  .  In  this  special  case,  the  reduced  Eqs.  30  and  31  are 
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identical  with  the  deep  water  relations  previously  obtained  by 
Manton  (15). 

The  dimensionless  significant  wave  height  follows  readily 
from  Eq.  20  as 

H  =  2  C2  (32) 

Hence,  by  specifying  functional  form  of  w(x,t),  various  special 
cases  can  be  examined  and  compared  with  other  available  para¬ 
metric  forecasting  theories  and  observational  data.  In  parti¬ 
cular, 


w(x,t)  =  H(x)  H(t)  (33) 

where  W(.)  =  Heaviside  unit  step  function,  corresponds  to  a 
duration-and  fetch-limited  generation  case  with  the  invariant 
solution 


C  +  In  (1-C)  =  ^  +  In  (1-^)  -  Y<)  (t  -  t<J) 


(34) 


x  =  x. 


*  'if'  £  (X,  T 


)  dT 


(35) 


where 


x  >  *o  ’  £  -  io  *  and»  Cq  =  £  (25c*  ic)  <  1  (36) 

Wind-wave  generation  under  conditions  of  uniform  wind  with 
unlimited  fetch,  i.e.,  w(t:)  *  W(t),  is  represented  by  the  solution 
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of  Eq.  30  with  d/dt  -  3/3t 


C  +  In  <1  -  C)  =  ^  +  In  (i  -  -  y0  (t  -  t^  (37) 

where 


t  >  to  ;  and,  ^  =  C  (t^)  <  1  (38) 

Finally,  the  solution  under  conditions  of  uniform  wind 
with  unlimited  duration,  i.e.,  w(x)  «  »(x),  also  follows  easily 
from  Eq.  30  with  d/dt  -  (C/2)  d/dx  as 

C  (1  +  h  C)  +  In  (1  -  C)  =  ^  (1  +  %  C0) 

+  In  (1-  C<))  -2yq  (x  -  x<,)  <39> 


where 


x  >  ajo  ;  and,  ~  £  (Xq)  <  1 

Noting  that  all  of  these  solutions  are  identical  to  Manton's 
(15),  we  can  proceed  now  to  compare  them  with  various  other 
forecasting  theories  and  observational  data  readily  available 
in  the  literature.  First,  the  fetch-limited  case  with  the  null 
boundary  condition  C-  ■  0,  i.e., 
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c  (1  +  %  C)  +  In  (1  -  C)  -  -2yo  X  (x  >  0) 


H  *  2  ^  C2 


where,  following  Manton  (15),  6*  1.31  x  10  2  and  Y0  =  2.5  x  10  , 

will  be  compared  with  the  parametric  relations  of  Hasselmann 
et  al.  (10): 


C  -  5.5  x  10"2  x0,3 


H  *  1.54  x  10~3  x0*5 


the  Sverdrup-Munk-Bretschneider  (SMB)  formulae  in  the  Shore 
Protection  Manual  (SPM)  (25) : 


C  =  1.2  tanh  (0.077  x0*25) 


H  =  0.283  tanh  (0.0125  x0,42) 


and  the  relations  by  Wilson  (30) 


C  *  5.7  x  10-3  x0,3 


H  =  2.5  x  10"3  x0,4 


The  wind-wave  parameters  represented  by  Eqs.  41  through  48 
are  shown  in  Figure  1  together  with  the  region  of  field  data 


obtained  from  Figure  15  of  Wilson  (30).  Evidently,  the  general 
character  of  the  results  generated  here  compares  favorably  with 
the  field  data  and  other  theories  or  empirical  curves  in  question. 
So  far  as  the  principal  purpose  of  any  parametric  approach  is 
concerned,  apparent  differences  observed  between  various  fore¬ 
casting  relations  can  be  considered  to  be  within  acceptable 
limits.  There  are,  however,  certain  points  of  concern  about 
each  forecasting  relation.  Specifically,  the  equations  of  the 
present  study  provide  estimates  with  inherent  upper  bounds  such 
that  C  <  1,  and  H  <  0.229.  These  bounds  appear  to  be  slightly 
underestimated  as  compared  to  some  field  data.  On  the  other 
hand,  parametric  relations  of  Hasselmann  et  al.  (10)  and  Wilson's 
(30)  empirical  curves  are  likely  to  overestimate  the  wind-wave 
parameters  at  large  fetches;  and, SMB  curves  overestimate  the 
same  parameters  at  small  fetches. 

The  duration-limited  case  with  the  null  initial  condition 
Cq  =  o,  i.e.. 


C  +  In  (1  -  C)  =  -yo  t  (t>0) 

(49) 

H  -  2  6^  C2 

(50) 

can  also  be  compared  with  the  relations  based  on  the  parametric 
approach  of  Hasselmann  et  al.  (10): 


C  *  9.17  x  10“3  t0,43 

(51) 

H  *  7.77  x  10"5  t5/7 

(52) 
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and  the  forecasting  curves  of  Sverdrup  and  Munk  (28),  which  are 
reproduced  here  in  Figure  2  together  with  the  region  of  data 
given  in  Figure  7  of  Ref.  28.  Also  shown  in  Figure  2  are  the 
wind-wave  predictions  based  on  Eqs.  49  through  52.  Again,  the 
overall  comparison  is  favorable.  Particular  differences  observed 
are  of  the  same  nature  as  those  just  mentioned  in  the  fetch- 
limited  case,  i.e.,  C  <  1,  H  0.229  for  the  predictions  of  the 
present  study,  which  are  likely  to  be  underestimates  for  long 
durations;  the  parametric  relations  of  Hasselmann  et  al.  (10) 
overpredict  C  and  H  tor  long  durations;  and,  the  Sverdrup  and 
Munk  theory  is  likely  to  overestimate  the  same  parameters  for 
short  durations. 

FORECASTING  IN  FINITE  WATER  DEPTH 

Uniform  Depth  Case.  -  Letting  =  0,  corresponding  to  a 
uniform  water  depth  d  =  const.,  and  with  0  *  @w,  being  coincident 
with  one  of  the  reference  axes,  the  governing  equations  take  the 
form 

dk/dt  =  (3/3t  +  V  3/3x)  k  *  (Sj  -  S2)  S^1  (53) 

dx/dt  =  V  =  n  C  (54) 

where  S^,  S2,  S^,  and  C  are  defined  as  before.  The  w^d-*  ve 
parameters  follow  from  Eqs.  19  and  20  as 

(55) 

(56) 
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T  «  {k  tanh  k  d)“* 


H  »  4 


<a/”n 


v“3  dv)^ 


We  can  now  proceed  to  consider  the  numerical  simulation 
of  a  set  of  three  wind-wave  generation  situations  corresponding, 
respectively,  to  fetch-,  duration-,  and  fetch-  and  duration- 
limited  wind  fields  with  invariant  direction.  In  all  these 
simulations,  we  take  Cf  =  0.01  as  a  characteristic  value  which 
allows  us  to  compare  the  results  generated,  in  particular,  under 
fetch-limited  conditions  with  those  presented  in  SPM  (25).  For 
the  remaining  cases  here  and  in  the  following  section  where  we 
will  examine  the  more  general  case  of  a  variable  bathymetry, 
there  is  no  parametric  results  available  for  comparison. 

The  fetch-limited  situation  corresponds  to  setting  d/djt  = 

V  d/dx  in  Eq.  53.  Hence,  the  equation  to  be  solved  becomes 

dk/dx  =  (Sj  -  S2)  /  (n  C  S4)  (57) 

where  w(x)  =  B  (x) ;  and,  T  and  H  are  given  by  Eqs.  55  and  56. 

In  this  case  and  in  all  subsequent  simulations,  null  boundary 

_]L 

and/or  initial  conditions  will  be  employed  so  that  ^  (=  k^) , 
T^,  and  are  all  zero  at  x  and/or  £  =  0.  Hence,  the  duration- 
limited  case  with  w(t)  *  B(£)  corresponds  to  the  solution  of 

dk/dt  =  (Sx  -S2)  S l1  (58) 

and  the  duration-and  fetch-limited  case  with  w(x,t)  *  W(x)  B(£) 
is  given  by  the  solution  of  the  full  set  of  Eqs.  53  and  54 
The  latter  can  also  be  expressed  in  terms  of  the  integral 
quadratures 


'•  .  /  5.  >v  r 


»**»•>'*  *■  1 
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(59) 


k  =  kD  +  f~  (S1  -  S2)  S'1  dr 


x  =  3^  +  f  n(k,x;d)  C(k,x;d)  dt 


(60) 


Many  sophisticated  numerical  algorithms  exist  for  the 
numerical  simulation  of  these  nonlinear  ordinary  differential 
equations.  However,  the  results  to  be  discussed  in  the  following 
were  generated  with  the  simplest  possible  scheme  by  approximating 
the  left-hand  sides  by  forward  finite  differences,  and  the  right- 
hand  sides  with  the  values  of  variables  k,  x,  £,  evaluated  at 
a  preceding  step.  Illustrating  this  with  Eq.  (58)  as  an  example, 
we  have  then 

ii.i  ■  h  *[<si  ‘  s2>  sZl]i4i  <6l> 

Since  the  null  initial  and/or  boundary  conditions  imply  that 
-+  oo  (or,  C_  *  0) ,  each  case  is  initiated  with  the  corresponding 
closed-form  solution  in  deep  water,  i.e.,  Eqs.  34  through  40 
whichever  is  applicable,  and  the  computation  via  closed- form 
solutions  are  continued  until  the  parameter  k  d  becomes  small 
enough  for  the  effects  of  transitional  water  depths  to  be  accounted 
for.  In  all  cases  here,  the  transitional  criterion  used  is 
k  d  =  10ff,  after  which  the  solution  is  continued  numerically.  It 
is  noted  that  for  the  duration-  and  fetch-limited  case  described 
by  Eqs.  53  and  54  or,  Eqs.  59  and  60,  it  is  expected  that  the 
local  field  is  initially  duration-limited  and  hence  grows  until 
it  becomes  fetch-limited.  The  boundary  between  the  fetch-limited 
region  and  the  duration-limited  region  is  given  by  Eqs.  59  and  60 


18 


■ 111  1 1 1  1 


with  x  =  t©  =  0  or,  equivalently,  from  a  comparison  of  the 
solutions  corresponding  to  the  fetch-1  imi  ted  and  duration- 
limited  cases.  In  other  words,  once  the  solutions  under  fetch- 
limited  ard  duration-limited  conditions  are  obtained,  the  fetch- 
and  duration-limited  case  follows  from  the  combination  of  the 
two  separate  cases  by  patching  them  together  in  a  straightforward 
manner. 

The  computations  for  the  fetch-limited  case  for  various 
d  values  are  shown  in  Figure  3  together  with  the  corresponding 
results  obtained  from  Eqs.  3-25  and  3-26  of  SPM  (25).  The 
general  agreement  between  the  two  sets  of  predictions  is  favor¬ 
able  except  for  short  fetches  where,  just  as  mentioned  previously, 
the  SPM  curves  tend  to  overestimate  the  parameters  T  and  H. 

Illustrated  in  the  subsequent  Figures  4  and  5  are  the  computations 
for  the  duration-limited  case  for  various  d  values  and  the 
fetch-and  duration-limited  case  for  d  =  0.3,  respectively.  There 
is  no  readily  available  observational  data  or  parametric  theories 
with  which  these  latter  results  could  be  compared. 

Variable  Depth  Case.  -  In  order  to  illustrate  the  application 
of  the  proposed  scheme  to  a  more  general  situation,  consider 
an  artifical  example  of  wind-wave  generation  over  a  circular 
basin  with  variable  depth 

d  =  1  +  10~2  (x2  +  £2)*  (x2  ♦  x2  1  95>  (62) 

Under  conditions  of  unlimited  wind-duration  with  w(x)  “  U(x) 
and  invariant  wind  direction  ©w  *  0  with  respect  to  +x  -  axis, 
the  governing  equations  follow  from  Eqs.  21,27,  and  29  as 

;jj 
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(63) 


dk/dt  *  J  .  £  k  -  (s1  -  S2  -  S3)  S'1 


d6/dt  =  -  (o/sinh  2  k  d)  m  .  V  d 

A 

(64) 

dx/dt  =  2  =  n  £  k 

(65) 

;  S^,  S2»  are  defined  by  Eqs.  23,  24,  26, 

i  as  in  Eq.  25  with  d(d)/dt:  =  0  ?  d;  and 

respectively; 

k  =  (cos  0,  sin  0) 

(66) 

A 

m  =  (-sin  0,  cos  0) 

(67) 

a  =  T  1  =  (k  tanh  k  d) ^ 

(68) 

The  parameters  T  and  H  are  given  by  Eqs.  19  and  20  as  before,  and 
6  is  with  respect  to  the  +x  -  axis. 

The  solution  in  this  case  proceeds  by  integrating  Eqs.  63 
and  64  along  the  rays  described  by  (65).  Hence,  using  forward 
differences  and  setting  As  =  nC  A£  as  a  stepsize  along  a  ray, 
the  computation  is  performed  in  the  following  order: 

0£+^  =0£-(om.£d/nC  sinh  2  k  d).  As  (69) 

x  =  x^  +  %  (cos  0£  +  cos  As  (70) 

i£+1  =  ££  +  %  (sin  0^^  +  sin  0£+1)  As  (71) 
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iii+1  s  hi  +  ^  “  ®3^i  + 

(st  ■"  ®2  ”  s3^i+i  (n  £  ^4^i+l^  ^s  (72) 

The  boundary  conditions  consist  of  the  null  values,  i.e.,  6  =  H  * 
T  *  0  for  any  ray  originating  along  the  semicircular  arc 
x  +  ^  =  95  (x  £  0)  corresponding  to  the  d  =  0.05  depth  contour. 

As  in  the  previous  uniform  depth  situations,  the  computation 
is  initiated  through  the  closed-form  solution  under  fetch-limited 
conditions,  i.e.,  Eqs.  41  and  42,  until  k  d  =  10ir,  after  which 
the  solution  proceeds  numerically  via  Eqs.  69-72.  The  results 
obtained  in  the  manner  just  described  with  As  »  0.2  are  presented 
in  Figure  6  through  8,  In  particular,  Figure  6  illustrates  the 
geometry  of  the  basin,  depth  contours,  wind  direction,  and 
computed  wave  rays,  which  remain  indifferent  to  refraction  over 
most  of  the  basin  except  towards  the  fetch-end  shoreline.  The 
contours  of  constant  T  and  H  are  shown  in  Figures  7  and  8, 
respectively. 


SUMMARY  AND  CONCLUDING  REMARKS 

Prediction  of  nearshore  wave  characteristics  generated  by 
variable  winds  was  considered  by  using  a  parametric  approach 
based  on  a  contracted  form  of  the  radiative  transport  equation. 

It  was  assumed  that  the  wind-disturbed  sea  surface  can  be  des¬ 
cribed  in  terms  of  a  depth-dependent  equilibrium  range  spectrum 
with  a  well-defined  low  wave  number  cut-off.  Among  all  possible 
processes  at  work  in  wind-wave  generation,  propagation,  and 
dissipation,  only  the  generative  effect  of  winds  and  the  dissipative 
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influence  of  turbulent  bottom  friction  were  taken  into  account, 
and  governing  equations  were  derived  to  describe  the  develop¬ 
ment  of  a  wind-wave  field  in  terms  of  the  integral  properties 
of  its  surface  spectrum  such  as  significant  wave  height  and 
characteristic  period.  The  governing  equations  were  solved 
first  for  various  simple  cases  under  deep  water  conditions, 
obtaining  results  which  compare  very  favorably  with  other  para¬ 
metric  theories  and  observational  data  available.  Consequently, 
the  solutions  were  numerically  extented  first  to  the  case  of 
wind-waves  in  uniform  water  depth  and  then  to  the  more  general 
case  of  a  circular  basin  with  variable  bathymetry.  Except  for 
the  case  of  wind-waves  generated  under  fetch-limited  conditions 
in  uniform  water  depths,  where  the  solutions  of  the  present 
study  compare  well  with  those  referenced  in  SPM  (25) ,  there  is  no 
other  corresponding  theory  or  readily  available  field  data  to 
verify  the  validity  of  the  remaining  set  of  numerical  results 
presented  in  Figures  4  through  8.  Hence,  these  results  are  of 
a  theoretical  nature  and  remain  to  be  substantiated  with 
observational  data  or,  by  generating  the  corresponding  integral 
properties  through  the  application  of  a  proven  discrete  spectral 
model  to  the  same  cases. 

Finally,  it  is  emphasized  that  the  proposed  parametric 
approach  is  an  approximate  one  involving  various  inherent  simpli¬ 
fying  assumptions  as  to  processes  which  operate  in  wind-wave 
generation.  Therefore,  the  proposed  approach  is  intended  for 
operational  uses  and  engineering  design  purposes  to  forecast  or 
hindcast  the  wind-wave  climatology  in  a  nearshore  zone  for  short- 
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or  long-term  applications  or,  to  fill  gaps  in  wave  data  base 
for  extreme  value  statistics.  Its  application  to  complex 
situations  such  as  hurricanes  involving  wind  fields  with  rapidly 
varying  directions  and  magnitudes  should  be  considered  only  with 
extreme  caution. 


Unlimited  Duration  in 


Comparison  of  Observational  Data  of  Wind-Wave 
Parameters  with  Various  Forecasting  Theories  for 
Conditions  of  Unlimited  Fetch  in  Deep  Water. 


Forecasting  of  Wind-Wave  Parameters  for  Conditions  of 
Unlimited  Duration  in  Uniform  Water  Depths:  a 
Comparison  of  SPM  Theory  with  the  Present  Study. 


Figure  4.  Forecasting  of  Wind-Wave  Parameters  for  Conditions 
of  Unlimited  Fetch  in  Uniform  Water  Depths  in 
Uniform  Water  Depths  Based  on  the  Present  Study. 


Figure  5.  Forecasting  of  Duration  and  Fetch-Limited  Wind-Wave 
Parameters  in  Uniform  Water  Depth  (d  *  0.3)  Based 
on  the  Present  Study.  “ 


Figure  7.  Contours  of  Scaled  Wave  Period  I. 


Figure  8.  Contours  of  Scaled  Significant  Wave  Height  H. 
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